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Global stability for an inverse problem in 
soil-structure interaction * 
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Abstract 

We consider the inverse problem of determining the Winkler subgrade 
reaction coefficient of a slab foundation modelled as a thin elastic plate 
clamped at the boundary. The plate is loaded by a concentrated force 
and its transversal deflection is measured at the interior points. We prove 
a global Holder stability estimate under (mild) regularity assumptions on 
the unknown coefficient. 
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1 Introduction 

The soil-structure interaction is an important issue in structural building design. 
The determination of the contact actions exchanged between foundation and 
soil is commonly approached by using simplified models of interaction. Among 
these, the model introduced by Winkler in the second half of the nineteenth 
century is one of the most popular in engineering applications |W) . In Winkler’s 
model, the foundation rests on a bed of linearly elastic springs of stiffness fc, 
fc > 0, acting along the vertical direction only. The springs are independent of 
each other, that is, the deflection of every spring is not influenced by the other 
adjacent springs. The accuracy of this model of interaction depends strongly 
on the values assigned to the subgrade reaction coefficient k. Ranges of average 
values of k are available in literature from extensive series of in-situ experiments 
performed on various soil types, ra, but these values are quite scattered 
and, in addition, they may vary significantly from point to point in the case of 
large foundations. Estimate of the coefficient k becomes even more difficult for 
existing buildings, since the soil on which the foundation is resting is not directly 
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accessible for experiments. For the reasons stated above, the development of a 
method for the determination of k is an inverse problem of current interest in 
practice. 

In this paper we consider the stability issue in determining the Winkler’s 
subgrade coefficient of a slab foundation from the measurement of the deflection 
induced at interior points by a given load condition. The mechanical model is 
as follows. The slab foundation is described as a thin elastic plate with uniform 
thickness h and middle surface coinciding with a bounded Lipschitz domain 
n C The plate is assumed to be clamped at the boundary 917, a condition 
that occurs when the slab foundation is anchored to sufficiently rigid vertical 
walls. A concentrated force of intensity / is supposed to act at an internal 
point Pq G Q. This load condition has the merit of being easy to implement in 
practice. According to the Winkler model of soil and working in the framework 
of the Kirchhoff-Love theory of plates, the transversal displacement w of the 
plate satisfies the fourth order Dirichlet boundary value problem 


div (div (^CV^ic))-1-A:u> =/i5(Po), in fl, 

(l.I) 

w = 0, 

on 9f7, 

(1.2) 

dw ry 

dn ~ 

on 9f7, 

(1.3) 


where C is the elasticity tensor of the material and n is the unit outer normal to 
917. Given the concentrated force f6{Po) and the coefficient k, k G L°°(f7), for 
a strongly convex tensor C G L°°(f7) the problem (ll.ll) - (ll.3l) admits a unique 
solution w G Hq{V,). 

The inverse problem in which we are interested in consists in studying the 
stability of the determination of the unknown subgrade coefficient k in (HID” 
(IE3D from a single measurement of w inside 17. It should be noted that the 
measurement of the transversal deflection at interior points of the plate can be 
easily carried out by means of no-contact techniques based on radar methodol¬ 
ogy (Ei]). 

Our main result states that, for C G IT^’°°(f7)ni/^+®(f7), for some 0 < s < 1, 
satisfying a suitable structural condition (see (13.31) 1. if Wi G i?o (^) solution 
to (ll.ll) - (|1.3l) for Winkler coefficient k = ki G D iJ'’(r7), 7 = 1,2, and if, 

for a given e > 0, 


||wi - W2\\L^n) < e/, 

(1.4) 

then, for every u > 0, we have 


\\ki-k2\\L^in^)<Ce^, 

(1.5) 


where 17^ = {x G 17 | dist{x, 9f7) > a} and the constants C > 0, /3 G (0,1) only 
depend on the a priori data and on a. 

It should be noted that one difficulty of the problem comes from the fact 
that the displacement w may change sign and vanish somewhere inside f7. See 
for instance the examples in [Gi] . |K-K-M| . fS^ . Therefore it is necessary to 
keep under control the possible vanishing rate of w. Thus, the key ingredients of 
the proof are quantitative versions of the unique continuation principle for the 
solutions to the equation div (div (^CV^w)) + kw = 0, precisely an estimate of 
continuation from an open subset to all of the domain (Propositionsand the 
Ap property (Proposition 13.5]) . Another useful tool is a pointwise lower bound 
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in a neighborhood of the point Pq where the force is acting for solutions to (ED 
fLemma l3.3|) . 

Let us mention that this method, essentially based on quantitative estimates 
of unique continuation, has some similarities, although with a different under¬ 
lying equation and with a different kind of data, with the one used in |A1] . for 
another inverse problem with interior measurements arising in hybrid imaging. 

The paper is organized as follows. Section [5] contains the notation, the 
formulation of the direct problem and a regularity result in fractional Sobolev 
spaces iProDOsition 12.311 . Section |2] is devoted to the formulation and analysis 
of the inverse problem. 

2 The direct problem 

2.1 Notation 

We shall denote by Br{P) the open disc in of radius r and center P. 

For any U CM.'^ and for any r > 0, we denote 

Ur = {x GU \ dist{x,dU) > r}. (2.1) 

Definition 2.1. regularity) Let LI be a bounded domain in Given k, a, 

with fc = 0,1,2,..., 0 < a < 1, we say that a portion S of dfl is of class 
with constants po, Mq > 0, if, for any P € S, there exists a rigid transformation 
of coordinates under which we have P = O and 

n ri Bp,,{0) = {x = {xi,X 2 ) G Bpg{0) I X 2 >'i){xi)}, 

where if is a (7^’“ function defined in Ip^ = {—po,po) satisfying 

m = 0 , 

'0'(O) = 0, when k>l, 
ll'0llc'‘.“(/po) < Mopo- 

When k = 0, a = 1, we also say that S is of Lipschitz class with constants po, 
Mo. 

We use the convention to normalize all norms in such a way that their terms 
are dimensionally homogeneous with the argument of the norm and coincide 
with the standard definition when the dimensional parameter equals one. For 
instance, the norm appearing above is meant as follows 

k 

2=0 


where 




a,Ip, 


(xi) - {yi)\ 

xi.yieip, \xi-yiU 

xWyi 


and denotes the i-order derivative of i). 
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Similarly, given a function m : fl i—!> K, where dfl satisfies Definition 12.11 and 
denoting by the vector which components are the derivatives of order i of 
the function u, we denote 


4L^{n)=Po { 


= Po' , A = 0,1,2,... 


Moreover, for fc = 0,1, 2,..., and s € (0,1), we denote 

IIwII//'=+«(n) = l|w||ff*^(n) +Po M^]s> 

where the semi-norm [ • ] s is given by 

\u{x)-u{y)\^ 


Q, 


- y\ 


‘2-\-‘2s 


dx dy 


( 2 . 2 ) 


For every 2x2 matrices A, B and every L G £(M^ x M^), 
following notation; 

{^A')ij — Lijf^iAf^i^ 

A ■ B = Aij Bij , 

|^| = (A.yl)i 

Finally, we denote by A'^ the transpose of the matrix A. 


we use the 

(2.3) 

(2.4) 

(2.5) 


2.2 Formulation of the direct problem 

Let fl C be a bounded domain whose boundary is of Lipschitz class with 
constants po, Mq and assume that 


|D| < MipI (2.6) 

We consider a thin plate Q x [—f,-!] with middle surface represented by 
and whose thickness h is much smaller than the characteristic dimension of D, 
that is h « pq. The plate is made by linearly elastic material with elasticity 
tensor C(-) G L°“(D, £(M^, M^)) with cartesian components CapjS satisfying 
the symmetry conditions 

CA = {CAf, (2.7) 

CA-B = A-CB, (2.8) 

for every 2x2 matrices A, B, and the strong convexity condition 

Coldp <Cd-d<6|d|2, (2.9) 

for every 2x2 symmetric matrix d, where positive constants. 

The plate is resting on a Winkler soil with subgrade reaction coefficient 

k G L°°{n), k>0 a.e. in fl. (2-10) 
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The boundary dQ is clamped and we assume that a concentrated force is act¬ 
ing at a point Pq € ^ along a direction orthogonal to the middle surface 
According to the Kirchhoff-Love theory of thin plates subject to infinitesimal 
deformation, the statical equilibrium of the plate is described by the following 
Dirichlet boundary value problem 


div (div (PV^zc)) + kw = 

in n, 

(2.11) 

w = 0, 

on 911, 

(2.12) 

dw n 

dn ~ 

on dVL, 

(2.13) 


where the plate tensor P is given by 

P=^C; (2.14) 

the subgrade reaction coefficient k satisfies 

fc 

0 < k(x) < a.e. in fl, (2-15) 

Po 

for some positive constant fc; the concentrated force is positive, i.e., 

/gR, />0; (2.16) 

w = w{x) is the transversal displacement at the point x € and n is the unit 
outer normal to cAl. 

We notice that the presence in (12.1111 and (12.1511 of the parameter po (which 
has the dimension of a length) allows for a scaling-invariant formulation of the 
plate equation. 

Proposition 2.2. Under the above assumptions, there exists a unique weak 
solution w G to (12.1111 - (12.1311 . which satisfies 

IkllH^(a) < Cf, (2.17) 

where the constant C > 0 only depends on h, Mq, Mi, ^q. 

Proof. The weak formulation of the problem (I2.11|I - (I2.13I1 consists in finding 
w G 77g(n) such that 

f PV^w • -I- f kwv = -^v{Po), for every z; G i7o(^)- (2-18) 

Jn JQ Po 

Let us notice that 

H^{n) C for every a < 1 (2.19) 

and, therefore, the linear functional 

F : R 

F{v) = ^v{Po) 

Po 
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is bounded and the symmetric bilinear form B(u, v) = PV^w ■ kwv 

is bounded and coercive on i?g(n) x iJQ(n). By Riesz representation theorem 
a solution to (12.181) exists and is unique. By choosing z; = tc in (12.181) . by (12.91) 
and using Poincare inequality, we have 

> f PV^w • V^w > (2.20) 

Po JO. Po 

where the constant C > 0 only depends on h, Mq, Mi, By (12.201) and the 
embedding (I2.19L the thesis follows. □ 


In the analysis of the inverse problem, we shall need the following regularity 
result when the coefficients of the plate operator belong to a fractional Sobolev 
space. 

Proposition 2.3 (i? ^-regularity). LetVL he a bounded domain with bound¬ 
ary of Lipschitz class with constants po, Mq, satisfying (12.61) . Given g € 
let w G be a solution to 

div {div {P\/‘^w)) = g, in Pi, (2-21) 

where P is given by (12.141) . with C satisfying (12.71) - (12.91) and, for some s € (0,1), 

l|C||w2.~(n) < -^2, (2.22) 

||C||jf2+s(Q) < Mq. (2.23) 

Then, for every a > 0, we have 

II'w|U4+«(q_^^j < C -I- Po\\g\\H-(o)'j , (2.24) 

where the constant C > 0 only depends on h, Mq, Mi, M 2 , Mq, s, a. 

Proof. When C gC°° {Pt) the estimate [)2.24p is a form of the well-known classi¬ 
cal Garding’s inequality. Under the less restrictive condition (12.221) . (12.231) . the 
proof of (12.241) can be carried out following the same path traced in the classical 
case i\M^ E) taking care to control the lower order terms by means of M 2 
and Mq. We omit the details. □ 


3 The inverse problem 

In order to derive our stability result for the inverse problem we need further a 
priori information. 

Concerning the point Pq of the plate in which the concentrated force is 
acting, we assume that 

dist{PQ,dP,) > dpQ, (3.1) 

for some positive constant d. On the elasticity tensor C = {Cap-ys} we further 
assume the stronger regularity (12.221) . (12.231) and, moreover, we introduce a 
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structural condition. Precisely, denoting by cq = Cim, oi = 4 Ciii 2 , 0,2 
2 C 1122 + 4(71212, 0,3 = 4(72212, 0,4 = (72222 and by S{x) the following matrix 


/ 

Oo 

Ol 

02 

03 

04 

0 

0 

\ 


0 

ao 

Ol 

02 

03 

04 

0 



0 

0 

Oo 

Ol 

02 

03 

04 



4ao 

3ai 

202 

03 

0 

0 

0 



0 

4ao 

3ai 

2a2 

03 

0 

0 



0 

0 

4oo 

CO 

202 

03 

0 


V 

0 

0 

0 

4ao 

3ai 

202 

03 

/ 


we assume that 

^{x) = 0, for every x G ft, 

where 

'D{x) = —|detS'(x)|. 
oq 


(3.2) 


(3.3) 

(3.4) 


Let us recall that condition (13.3|) includes the class of orthotropic materials and, 
in particular, the isotropic Lame case, see |M-R-V) . Concerning the subgrade 
reaction coefficient k, we require the additional regularity 


Pq < 4 - ( 3 - 5 ) 

Po 

Remark 3.1. Let us emphasize that the assumption k € is not merely 

a mathematical technicality, but it can be grounded on realistic mechanical 
considerations. If, for instance, k is piecewise constant and is represented as 


k{x) = ^ kjXEj {x), for every x € (3.6) 

1=1 

where kj G R and Ei,...,Ej is a partition of into disjoint subsets of finite 
perimeter (in the sense of Caccioppoli, that is XEj G i3fo(R^) for every j), then 
k belongs to i7®(]R^) for every s, 0 < s < ^. In fact one has 

[krs<Cs\\k\\f^([ ' 7 / mT, 

\7r2 / VJb2 / 

for every k G L^(]R^)nL°°(R^)(7iIfo(R^). Here Cs only depends on s G (0, |) and 
Jg 2 \Dk\ denotes the total variation of k. For a proof see |M-P1 formula (2.15)] 
and also [Gil Remark 1.16] for the convergence properties of the mollifications 
of BV functions. 

In particular, if k is of the form (13.61) and we assume 

= f \PXE, I < Vpo, for every j = 1,..., J, 

for a given R > 0, then we obtain 
[fc]2 < 

Hereinafter, we shall refer to h, d, Mq, Mi, M 2 , M 3 , ^ 0 , k, s as the a priori 
data. 
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Theorem 3.2. Let fl be a bounded domain in with boundary of Lipschitz 
class with constants po, Mq, satisfying (|2.6I) . Let P given by (12.141) . with C G 
n satisfying (12.71) - (12.91) . (12.221) . (12.231) for some s G (0,1), 


and ( 13 . 31 ) . Let Pq G Q satisfying ( 13 . 11 ) . 

Given / > 0, let Wi G i = 1,2, be the solution to 

{ div (div (FV'^Wi)) + kiWi = in 11, (3-7) 

Po 

Wi = 0 , on dLl, ( 3 - 8 ) 

= 0 , on dn, ( 3 . 9 ) 

for ki G L°°{n) n H^{n) satisfying (| 2 . 15 l) and ( 13 . 51 ) . 

If, for some e > 0, 

\\wi-W 2 \\L^(a) <ef, (3.10) 

then for every a > 0 we have 

\\ki - fc 2 ||L 2 (n^p„) < —( 3 . 11 ) 


where the constants C > 0 and ft G (0,1) only depend on the a priori data and 
on a. 

As is obvious, the above stability result also implies uniqueness. Indeed, by 
the following arguments it is easily seen that, under the above stated structural 
conditions on C (iT^-dSa. uniqueness continues to hold by merely assuming 
fc G and C G 

Let us premise to the proof of Theorem 13.21 some auxiliary propositions 
concerning quantitative versions of the unique continuation principle (Lemma 
and Propositions 13.41 and 13.51 belowL 

Lemma 3.3. Let LI be a bounded domain with boundary dLl of Lipschitz class 
with constants po, Mq, satisfying (12.61) . Let Pq G LI satisfying p.ll) . Let P 
given by (j2.14p . with C satisfying (12.71) - (12.9|) . and let k and f satisfy (I2.15|) . 
(I2.16p . respectively. Let w G Hq{ 1}) be the solution to (12.111) - (12.13p . There 
exists a > 0, only depending on h, d, Mq, Mi, Ci; such that 

w{x)>Cd^f, yx € B 2 apoiPo), (3.12) 

where C > 0 only depends on h, Mq, Mi, 

/ > C'cr^d^pQ||w;||^ 2 (f 2 ), for every a, 0 < cr < u (3.13) 

i^o)\^fTpQ {Po) 

where C > Q only depends on h, Mq, Mi, Cij k. 

Proof. By (12.91) and (I2.18p . we have that for every v G Hq{LI) 

f\v{Po)\ < C'll'a||ff 2 (n)||w||// 2 (f 2 ), (3.14) 

so that 

ll^(7o)||ff-2(n) = sup |'y(Pb)| < ^Ikllff 2 (n), (3.15) 
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where C > 0 only depends on h, k. Since Bdpo{Po) C n by (13.11) . we have 

ll<^(-Po)||ff-2(n) > ll<^(^o)||ir-2(Brfp^(Po)) > Cd, (3.16) 

where C > 0 is an absolute constant. From (I3.15L (13.1611 it follows that 

Iklk^(n) > Cdf, (3.17) 

where C > 0 only depends on h, k. By (12.911 . (12.1811 and Poincare inequality, 
we have 

w{Po) > y|k||^2(n), (3.18) 

where C > 0 only depends on h, Mq, Mi, By (I3.17|) . p.l8l) and by the 
embedding inequality (12.191) . we have 

w{Po) > Cd\\w\\H2(n), (3.19) 

w{Po) > cod||w||co,c(Q), (3.20) 

where C > 0 and cq > 0 only depend on h, Mq, Mi, ^q, Ci, k. Let 

" = (3.21) 

Let us notice that, by this choice of u, dist{Po, dV,) > iapo and, recalling (I3.2()L 
we have 

w{x) > w{Po) - |w(a:) - w(Po)| > w{Po) - (2CT)“||u;||c;o,c.(n) > 

for every x G B2apoiPo)- (3.22) 

Choosing a = ^, (13.121) follows from (13.171) . (13.1811 and (13.221) whereas (13.131) 
follows, restricting to the annulus B^apoiPo) \ BrrpoiPo), which is contained in 
B 2 Wpo{Po) for a < a, from (13.1911 and (13.2211 . □ 


Proposition 3.4 (Lipschitz propagation of smallness). Let U be a bounded 
Lipschitz domain o/R^ with constants po, Mq and satisfying |17| < Mip^. Let 
w G H'^{U) be a solution to 

div{div{FV'^w)) + kw = 0, in U, (3.23) 

where P, defined in (12.141) . satisfies (12.711 - ()2.9I) and (I2.22p in Lf, and k satisfies 
(I2.15|l in Lf. Assume 


\\w\\l^U) 

There exists a constant ci > 1, only depending on h, M 2 , fo and k, such that, 
for every t > 0 and for every x G Lfc^Tpo! we have 



> Ct 


where Ct > 0 only depends on h, Mq, Mi, M 2 , k, r and on N. 


(3.24) 
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The proof of the above proposition is based on the three spheres inequality 
obtained in c 'WW\ . 

Proposition 3.5 (Ap property). In the same hypotheses of Proposition 
there exists a constant C 2 > 1, only depending on h, Mq, Mi, M 2 , k, such 
that, for every r > 0 and for every x S Uc^rpo, we have 


1 

\ Btpo (3^)1 



1 

\Btpo{x)\ 



<B, 


(3.25) 


where B > 0 and p > 1 only depend on h, Mq, Mi, M 2 , k, t and on N. 


The proof of the above proposition follows from the doubling inequality 
obtained in |dC-L-M-R-V-W| . bv applvine the arguments in 

Proof of Theorem \3.‘A If e > 1, then the proof of (13.111) is trivial in view of 
(12.1511 . Therefore we restrict the analysis to the case 0 < e < 1. 

The difference 

w = wi — W 2 (3.26) 

of the solutions to (I3.7I) - (I3.9I) for i = 1,2 satisfies the boundary value problem 

( div (div (FV^w)) + k 2 W = (k 2 — ki)wi, in 11, (3.27) 

J re = 0, on on, (3.28) 

[ 1 ^ = 0 , on 511. (3.29) 

Obviously, it is not restrictive to assume that a < a, where a has been defined 
in (13.211) with a = 5 . We have 



kifwl<2{h+l2), 


where 

h= [ klw'^, 

•tOo-po 

l2= [ (div(div(PV2^^;)))^ 

By (|2.15ll and (I3.1QL we have 


(3.30) 


(3.31) 

(3.32) 


Ii < —< 

Po 


(3.33) 


By (|2.22L we have 




12 p« 


(3.34) 


with C > 0 an absolute constant. Let g = (^2 — ki)wi — k 2 W. Note that, by 

( 1 ^ 1 ^ . dUZl), (EUl), and dim . 


Po 


kf 


(3.35) 
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where C > 0 only depends on h, Mq, Mi, By applying Proposition 12.31 we 
have 

< Cfk, (3.36) 

with C > 0 only depending on h, Mq, Mi, M 2 , M 3 , ^ 0 , s, a. From the well- 
known interpolation inequality 


and recalling (13.361) and (13.1011 . we obtain 

Ik||ff4(n ) < Cfe^, 


(3.37) 


(3.38) 


where C > 0 only depends on h, Mq, Mi, M 2 , M 3 , ^oi k, s, a. From (I3.30L 
(I3.33p . (I3.34|) . (13.381) . it follows that 


/ (^2 - kifwl < , 

Po 


(3.39) 


where C > 0 only depends on h, Mq, Mi, M 2 , M 3 , ^O) k, s, a. 

Let us first estimate 1^2 — fci| in a disc centered at Pq- Notice that, by the 
choice of CT, flo-po D i? 2 CTpo(^b)j for every a < a. By applying (13.121) for w = wi, 
and by (13.3911 with a — a, we obtain 


[ (fc 2 -fci )2 < 

dB2^po(Po) ^ 0 “ 

where C > 0 only depends on h. Mg, Mi, M2, M3, ^ 0 , k, s and d. 
Now, let us control 1^2 — fci| in 

f^crpo = f^o-po \ B2apo ■ 


(3.40) 


(3.41) 


This estimate is more involved and requires arguments of unique continuation, 
precisely the ^p-property and the Lipschitz propagation of smallness. 

By applying Holder inequality and (13.3911 . we can write, for every p > 1, 


f (* 2 -fcl)^=/ \wi\i{k 2 - kif\wi\ P< 

*^Qcrpo 

- i L {k2-kifwl\ ([ {k2 - ki^lwil' 

\«/Qcron / \*^^cron 


< 


C 2 2s 

< —— J P gp(4-f-s) 

Po 


{k 2 - ki)'^\wi\ p-i 


(3.42) 


where C > 0 only depends on h, Mq, Mi, M 2 , M 3 , ^ 0 , k, s, a. 

Let us cover ftapo with internally non overlapping closed squares Qi (xj ) with 
center xj and side I = ^Tnaxf^jTcT }^ ~ f ’ '''^here ci and C 2 have been 
introduced in ProDOsition l3.4l and in Proposition l3.5l respectively. By the choice 
of I, denoting r = 


J .7 

f^o-po C u Ql{xj) C (J Bj.{Xj) C f^fpo \ BapoiPo), 

1=1 7=1 


(3.43) 
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so that 


f {k2-h)^\wi\ |wi| 

•^n^po Po •tn.^po Po 


kiT 


(3.44) 


By applying the ^p-property (13.251) and the Lipschitz propagation of smallness 
property ([323 to u> = u>i in 17 = \ B^pjPo), with r = ^ = 2 max{lci,e 2 } ’ 
and noticing that, for every j = 1,J, dist{xj^dU) > cir^ i = 1, 2, we have 


f Br{Xj) 


_ 2 




1 — 

(|bX)| t 


|Bp(a:j)| 


'n\B^,^(Po) 


K |2 


1 7 

P — 1 


(3.45) 


where i? > 0, p > 1 and Cr > 0 only depend on h, Mq, Mi, M 2 , k, a and 

ll“ill^i ^ ^ ^ 

on the frequency ratio T = -n—n- ——Such a bound can be achieved 

II’"iIIi,2{SJ\B_p^(P0)) 

as follows. Notice that, since fl\Bapg{Po) D i? 2 apo(a) \Bwpo{Po), by applying 
(13.1311 . we have 

P<7 (3.46) 

where C > 0 only depends on h, Mq, Mi, ^q, Cii k. By applying (13.1311 and 
(I3.17|l to estimate from below the denominator in the right hand side of (I3.45L 
by (12.61) and (13.441) . we obtain 


f (^2 - fcl)^|wi| P-i < 




C 


pg(d4/2)pli p6M4/2)pii 


(3.47) 


where C > 0 only depends on h, Mq, Mi, M 2 , k, d and a. By (13.421) and 
(13.4711 we have 

f ik2 - kif <, (3.48) 

where C > 0 only depends on h, Mq, Mi, M2, M3, ^O) k, s, d and a. Finally, 
by (13.401) and (13.481) . the thesis follows. □ 
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